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ABSTRACT 

NRAO 530 is an optically violent variable source and has been studied with multi- 
epoch multi-frequency high-resolution very long baseline interferometry (VLBI) ob- 
servations. NRAO 530 was monitored with the Very Long Baseline Array (VLBA) at 
three frequencies (22, 43 and 86 GHz) on 10 consecutive days in 2007 May during 
observations of the Galactic center (Sgr A*). Furthermore, analysis of archival data of 
NRAO 530 at 15 GHz over the last 10 years allows us to study its detailed jet kinemat- 
ics. We identified the compact component located at the southern-end of the jet as the 
VLBI core, consistent with previous studies. The 10-d monitoring data at the three 
high frequencies were shown to produce high-quality and self-consistent measurements 
of the component positions, from which we detected for the first time a two-dimensional 
frequency-dependent position shift. In addition, the repeated measurements also per- 
mit us to investigate the interday flux density and structure variability of NRAO 530. 
We find that it is more variable for the inner jet components than those further out. 
We obtained apparent velocities for eight jet components with /3 a pp ranging from 2 
to 26c. Accordingly, we estimated physical jet parameters with the minimum Lorentz 
factor of 14 and Doppler factors in the range of 14-28 (component f). The changes 
in the morphology of NRAO 530 were related to the motion of separate jet compo- 
nents with the most pronounced changes occurring in the regions close to the core. For 
NRAO 530, we estimated a position angle swing of 3? 4 per year for the entire inner jet 
(components d and e) . The non-ballistic motion and change of jet orientation make 
this source another prominent example of a helical and possibly "swinging" jet. 

Key words: galaxies: jets - quasars: individual: NRAO 530 - radio continuum: galax- 
ies. 



1 INTRODUCTION 

At a rcdshift of 0.902 (Junkkarinen 1984), NRAO 530 (also 
known as B1730-130) is a well-known blazar, which is clas- 
sified as an optically violent variable object. Strong and er- 
ratic broad-band variability has been observed in the radio 
(Bower et al. 1997), optical (Webb et al. 1988), and the 7-ray 
regime (Mukherjee et al. 1997). In 7-rays, it appeared as a 
relatively quiescent source after the launch of Fermi (Abdo 
et al. 2009), but started to flare from the end of 2010 Oc- 
tober (D'Ammando & Vandenbroucke 2010). Interestingly, 
Foschini et al. (2006) serendipitously detected a short hard 
X-ray flare on 2004 February 17 with the INTEGRAL satel- 
lite. 

On kpc scales, NRAO 530 exhibits two-sided lobes in 
the east-west direction. The western lobe is stronger than 
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the eastern one and is connected to the core with weak 
arch-structured emission (Hong et al. 2008). Centimetre very 
long baseline interferometry (VLBI) images showed a core- 
jet structure on pc scales with an oscillating trajectory ex- 
tending to the north from the assumed core (e.g. Shen et al. 

1997) . Space VLBI observations revealed brightness temper- 
atures of NRAO 530 significantly in excess of both, the in- 
verse Compton and the equipartition limit (Bower & Backer 

1998) . Superluminal motion of several jet components is de- 
tected with apparent velocities in the range of 10-40c (Bower 
et al. 1997; Jorstad et al. 2001; Feng et al. 2006; Hong et al. 
2008). 

Images of NRAO 530 obtained at 5 GHz in the early 
1990s indicated that the position angle (PA) of the jet com- 
ponents changed significantly (Hong et al. 1999). Our new 
images in 2007 reveal a jet structure for the inner region 
very different from those previously reported (e.g. in Feng 
et al. 2006). NRAO 530 has undergone two moderate flares 
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(Aller, Aller & Hughes 2009; Jenness, Robson & Stevens 
2010) since the dramatic flare around 1997 (Bower et al. 
1997; Feng et al. 2006). This stimulates us to also investi- 
gate the pc-scale jet kinematics and its possible relation to 
the long-term radio flux variations in NRAO 530. 

Following a brief explanation of the observations and 
data reduction in Section 2, we present in Section 3 the re- 
sults from Very Long Baseline Array (VLBA) observations 
at three frequencies (22, 43 and 86 GHz). The structure of 
the inner jet was found to have changed significantly when 
compared with previous VLBI studies. In order to better 
understand the jet kinematics, we therefore have made use 
of VLBA data at 2 cm (15 GHz), from the Monitoring Of 
Jets in Active galactic nuclei with VLBA Experiments (MO- 
JAVE) programme, which allows us to conduct a detailed 
study of structural variations with good time coverage. We 
briefly discuss in Section 4 the implications for our under- 
standing of jet morphology and kinematics. A summary of 
the results is presented in Section 5. 

Throughout this paper, we adopt the luminosity dis- 
tance to NRAO 530 £>l = 5.8 Gpc, 1 mas of angular sep- 
aration corresponding to 7.8 pc, and a proper motion of 
lmas yr _1 corresponding to a speed of /3 apP = 48.5 c (Ho 



NRAO 530 at 22.228 GHz in I 2007 May 15 



= 71kms~ 1 Mpc" 
2009). 



O.M = 0.27, r?A = 0.73, Komatsu et al. 



2 OBSERVATIONS AND DATA ANALYSIS 

During the global campaign on SgrA* in 2007 May, 
NRAO 530 was monitored as a calibrator with the VLBA 
from 2007 May 15 to 24 (2007.370-2007.395) at 22, 43 
and 86 GHz every day. NRAO 530 was observed as a fringe 
finder and was used to check the amplitude calibration of 
SgrA* (Lu et al. 2011). NRAO 530 was observed shortly 
before, during, and after the time of visibility of the tar- 
get source, SgrA*. Dual circular polarization at a recording 
rate of 512 Mbps [8 intermediate frequency (IF) channels, 
16 MHz per IF, and 2 bits per sample] was recorded at 
each station. The individual VLBI scans on NRAO 530 had 
a duration of 3min at 22 GHz and 4min at 43 and 86 GHz, 
respectively. Five scans at each frequency were observed, 
resulting in a total time on source of ~ 15 min at 22 GHz 
and ~ 20 min at 43 and 86 GHz, respectively. An example 
of the uv coverage at 22 GHz is shown in Fig. 1. The data 
analysis was performed within the AIPS software in the usual 
manner. Opacity corrections were made at all three frequen- 
cies by solving for receiver temperature and zenith opacity 
for each antenna. In addition, we have also used the VLBA 
observations at 15 GHz from the "2 cm Survey" and the fol- 
low up MO J AVE 1 program between 1999 and 2009, spanning 
10 years. These data were provided as automatically self- 
calibrated uv FITS file, which we re-do the mapping and 
self-calibration in difmap before model-fitting. 

In all cases of the model fitting, circular Gaussian com- 
ponents have been chosen in order to simplify and unify 
the comparison. The final fitting of the jet components was 
stopped when no significant improvement of the reduced 
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Figure 1. An example of uv coverage at 22 GHz on 2007 May 
15. 



Xv values was obtained. The formal errors of the fit pa- 
rameters were estimated by using the formalisms described 
in Fomalont (1999). The position errors were estimated by 
cr r = ,f "" s d , where u rms is the post-fit rms, and d and 7 D eak 

2 'peak 

are the size and the peak intensity of the component. In 
case of very compact components, this tends to underes- 
timate the measurement error. We therefore included an 
additional minimum error criterion according to the map 
grid size, which roughly matches one-fifth of the beam size 
at each frequency and corresponds to 0.06 mas at 15 and 
22 GHz, 0.03 mas at 43 GHz, and 0.02 mas at 86 GHz. For 
the weighted mean of parameters from the 10 consecutive 
epochs, we have chosen the max (^/l/E^y , -^7=) as the stan- 
dard error for the few cases when the former is smaller than 
the latter, with a being the standard deviation. By doing 
so, we have assumed stationarity of the corresponding pa- 
rameter on a daily time-scale. For the spectral analysis, we 
estimate absolute flux-density calibration errors of 3-5 per 
cent at 22 and 43 GHz and 20 per cent at 86 GHz. Our 
10-d observations have revealed a consistent jet structure 
at all three frequencies. Cross-identification of jet compo- 
nents is straightforward and unambiguous. The components 
seen at 15 GHz were cross-identified based on their flux den- 
sity evolution, core separation, PA, and full width at half- 
maximum (FWHM). We noted that although the archival 
data at 15 GHz have a good time resolution, large time gaps 
between 2000 and 2002 pose some difficulties in tracing the 
weak jet components (see Section 3.4). 



3 RESULTS 

In Fig. 2, we show as an example one of the resulting 
clean maps of NRAO 530 at 15 GHz at epoch 2007 June 
10 (2007.441) and at 22, 43, and 86 GHz at epoch 2007 May 
15 (2007.370). Maps for all epochs at all frequencies are in- 
cluded as Supporting Information in the online version of 
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Figure 2. Example of clean maps of NRAO 530 at (a) 15 GHz, (b) 22 GHz, (c) 43 GHz, and (d) 86 GHz. See Table 1 for the map 
parameters, and the Supporting Information section for all images. 



this paper (Figs 14 and 15). Table 1 lists the map parame- 
ters and Table 2 summarizes the epoch of observation, the 
identification of components, and the model fit parameters 
and their uncertainties (all errors are 1 a throughout) of 
the maps shown in Fig. 2 (only a small portion of each ta- 
ble is shown here; the full tables are available online in Ta- 
bles 8 and 9 - see Supporting Information). On 2007 May 15 
(2007.370), the sparse uv coverage at 86 GHz did not per- 
mit us to model the source structure by more than a single 



Gaussian component, and we will not consider this data set 
for further analysis. 

On VLBI scales, NRAO 530 is characterized by a one- 
sided jet following a curved trajectory of ~25mas to the 
north (e.g. Romney et al. 1984; Bondi et al. 1996). We find 
no naming convention for the VLBI components in the lit- 
erature that we could follow for NRAO 530. We therefore 
labelled the modelled components in an alphabetical order 
with increasing core separation, starting from the core (c) 
to the outermost jet component (j). A very faint jet feature, 
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Table 1. Description of VLBA images of NRAO 530 shown in Fig. 2. The columns list the epoch of observation, the observing frequency, 
the peak flux density, the parameters of the restoring elliptical Gaussian beam (the FWHM of the major and minor axes and the PA 
of the major axis), and the contour levels of the image, expressed in percentage of the peak intensity. See the Supporting Information 
section for the full table. 









Restoring Beam 




Epoch 


V 


"Speak 


Major 


Minor 


P.A. 


Contours 




(GHz) 


(Jy beam -1 ) 


(mas) 


(mas) 


(°) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2007.441 


15 


1.20 


1.34 


0.568 


-5.7 


-0.35, 0.35, 0.7, 89.6 


2007.370 


22 


1.32 


1.19 


0.344 


-12.8 


-0.3, 0.3, 0.6, 76.8 


2007.370 


13 


1.33 


0.537 


0.193 


-7.7 


-0.3, 0.3, 0.6, 76.8 


2007.370 


86 


0.74 


0.2 


0.2 





-1, 1, 2, 64 



Table 2. Parameters of the Gaussian model fit components of the images shown in Fig. 2. See the Supporting Information section for 
the full table. 



Epoch 


Id. 


Flux 
[mJy] 


Core Separation 
[mas] 


P.A. 

[degree] 


Size 
[mas] 






v = 


15 GHz 






2007.441 


c 


1138.4±60.1 


O.OOiO.OO 


O.OiO.O 


0.07±0.01 




d 


545.9±20.6 


1.17±0.06 


-8.6±2.9 


0.31±0.01 




e 


310.9±47.4 


1.52±0.06 


15.4±2.3 


0.62±0.09 




f 


109.6±19.5 


3.29±0.06 


27.0±1.0 


0.74±0.11 




g 


67.1±19.4 


5.55±0.35 


5.0±3.6 


2.49±0.70 




i 


139.5±22.7 


8.87±0.12 


13.5±0.8 


1.52±0.23 




j 


121.9±37.0 


24.19±0.65 


0.9±1.5 


4.33±1.30 






v = 


22 GHz 






2007.370 


c 


1332.9±23.6 


O.OOiO.OO 


O.OiO.O 


0.08±0.01 




d 


385.9±40.5 


1.10±0.06 


-9.5±3.1 


0.28±0.02 




e 


377.8±72.1 


1.52±0.08 


9.6±2.9 


0.87±0.15 




f 


99.1±26.5 


3.41±0.10 


26.1±1.6 


0.81±0.20 




g 


29.5±11.2 


4.52±0.22 


7.0±2.8 


0.32±0.09 




i 


138.3±44.2 


8.75±0.30 


13.4±2.0 


1.99±0.61 






v = 


43 GHz 






2007.370 


c 


1258.0±60.7 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




n 


272.9±22.4 


0.38±0.03 


-3.2±4.5 


0.08±0.01 




d 


221.2±27.7 


1.31±0.03 


-10.4±1.3 


0.21±0.02 




e 


250.3±40.2 


1.56±0.05 


7.5±2.0 


0.70±0.11 




f 


56.1±20.2 


3.44±0.12 


24.4±2.0 


0.68±0.24 






v = 


86 GHz 






2007.370 


c 


1214.4±186.3 


O.OOiO.OO 


0.0 ±0.0 


0.16±0.02 



which faded away after 2002 is labeled as x. An inner jet 
component, which was not resolved at frequencies 22 GHz 
in 2007 May, is named as n and the newly ejected compo- 
nents after 2008 is named as n' . In our analysis, we identify 
the southernmost component (component c) to be the sta- 
tionary core of the radio emission, because of its flat radio 
spectrum (see Section 3.1) and high compactness (see Ta- 
ble 2). This is in agreement with the identification scheme 
proposed by Jorstad et al. (2001), Feng et al. (2006), and 
Chen, Shen & Feng (2010). 

As can be seen in Fig. 2, our three-frequency data sets 
of 2007 revealed a very consistent jet morphology within 
the central ~ 10 mas, with richer structure seen at lower 
frequencies. The innermost jet component (n) detected at 
43 and 86 GHz was not seen at 22 GHz, which may be due 
to the insufficient angular resolution at this frequency and 
the relative weakness, when compared with adjacent com- 
ponents. The inner jet extends slightly to the north-west 
direction (PA ~ —10°) out to a core separation of ~ 1.5 mas 



and then bends sharply by ~ 90° towards a PA of 25° at a 
core separation of 3.5 mas. From here, the jet bends gently 
northward following a curved path with underlying diffuse 
emission (see Fig. 2). Finally, it fades away at a core sep- 
aration of ~ 9 mas to a PA ~ 10°. The jet component (i) 
located at ~25 mas north of the core is not seen at the three 
highest frequencies, but is clearly visible at 15 GHz across 
all the epochs (Table 2). This is probably due to the lower 
resolution at the 15 GHz and the steep spectrum of this jet 
component. 

3.1 Component spectrum 

In order to study the spectra of the VLBI components, we 
also used the 15-GHz data of 2007.441, taken close in time to 
our multifrequency observations in 2007 May. The spectra 
of the components are shown in Fig. 3 based on the flux 
densities obtained by model-fitting of Gaussian components 
to the visibility data. The spectrum of the core component 
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Figure 3. Component spectra of NRAO530 in 2007.4. Left: Spectrum of the core. The SSA fit is shown as a solid curve. Right: Spectra 
of the jet components. 



c, which is located at the south end of the VLBI jet, is shown 
in Fig. 3 (left). The presence of an unambiguous turnover 
and the flattest optically thin spectral index (—0.19 ± 0.03 
between 43 and 86 GHz) among all the components allow us 
to identify this component as the compact VLBI core. 

By assuming that the spectral turnover is due to syn- 
chrotron self-absorption (SSA; Jones, O'dell & Stein 1974), 
the convex spectrum of the core can be fitted by the follow- 
ing equation: 

S v = So^ 2 ' 5 [l - exp(-r s /, Q - 2 ' 5 )], (1) 

where v is the observing frequency in GHz, So is the intrin- 
sic flux density in Jy at 1 GHz, r s is the SSA opacities at 
1 GHz, and a is the optically thin spectral index. The best- 
fit parameters are a — —0.22 ±0.01, So = 2.2 ±0.1 mjy and 
r s = 1342.7 ± 14.4. The fitted spectra are shown in Fig. 3 
(left). 

The magnetic field B of a homogeneous synchrotron 
self-absorbed source component can be calculated via 

B syn = 10 -B 6(a)l/ « ax e± Sm l x &/{l + z), (2) 

where v mBX is the peak frequency in GHz, 6 the source an- 
gular size in mas, S ma ,x the peak flux density in Jy and 
b(ot) a tabulated function of the spectral index ct (Marscher 
1983). Assuming similar Doppler factors S for the core and 
jet component f (5 = 14.1, see Section 3.4), we obtain 
B syn = 76.1 mG with the following parameters: 6(a) = 1.8, 
Smax = 1.339 Jy, //max = 27.5 GHz, 9 = 0.09 mas. This is not 
in agreement with the magnetic field estimation in Feng et 
al. (2006), who obtained 8mG for the core component by 
assuming 5 = 1. 

The spectra of the outer jet components are shown 
in Fig. 3 (right), along with their spectral fitting. Com- 
ponent n is most likely self-absorbed at frequencies below 
43 GHz. From 43 and 86 GHz, the extrapolated flux density 
at 22 GHz is ~ 0.9 Jy. Without self-absorption, such a bright 
feature should have been seen at 22 GHz at a core distance 
of ~ 0.4 mas (see also Fig. 5). Other jet components (d, e, f, 
g and i) displayed steep spectra, except for some indication 
of absorption for component e at 15 GHz. Further observa- 
tions with wider frequency coverage for component e would 
be needed to clarify this. The corresponding optically thin 
spectral indices a (S v oc v a ) for these jet components are 
a n = -1.90±0.19, a d = -0.93±0.09, a c = -0.51±0.21, at 



= -0.71 ± 0.03, Q g = -1.89 ± 0.92, and an = -0.53 ± 0.51, 
where the index labels each component. 



3.2 Frequency-dependence of component positions 

The 10-d repeated measurements at 22, 43 and 86 GHz allow 
us to investigate the frequency-dependent position shift of 
the jet components. In Fig. 4, we show the positions of the 
components detected at any two of these three frequencies. 
By using the additional data at 15 GHz (epoch 2007.441), 
the spectra of the core and components d, e and f are shown 
as insets in Fig. 4. The time averaged component positions 
are summarized in Table 3. 

As it can be clearly seen in Fig. 4, the relative positions 
of three components ( d, e, and f) are systematically dis- 
placed between 22 and 43 GHz. These are components whose 
positions are most reliably measured in our experiments. For 
all three components, the core separation at 43 GHz is larger 
than that at 22 GHz, at a > 3 a level. For the components 
d, e and f, the position shift between 22 and 43 GHz is 0.17, 
0.12 and 0.21 mas, respectively (see Table 3). The measured 
position of component d at 86 GHz is consistent with the 
lower frequency data in the sense that its core distance at 
86 GHz is larger than that at lower frequencies. For compo- 
nent n, we are unable to compare the position shift between 
43 and 86 GHz, limited mainly by the large position scatter 
at 86 GHz. At the same time, the observations at 15 GHz 
also provide consistency checks for the components d, e and 
f. To further verify this position shift effect, we show in Fig. 5 
radial slices made along PA of —10° for the inner ~ 2 mas of 
the jet. One can clearly see that the peak of the component 
d is displaced between 22, 43 and 86 GHz. For the inner- 
most jet component n, however, it is impossible to reliably 
model it at frequencies ^ 22 GHz, and its position measured 
at 86 GHz suffers from large scatter. 

We also used spectral index mapping to pin down this 
effect. In Fig. 6, we show a spectral index map constructed 
between 22 and 43 GHz. Prior to making the spectral in- 
dex map, we have corrected the misalignments between the 
two maps by shifting one map relative to the other. Such 
an alignment is necessary as suggested by the appearance of 
obvious false features in spectral index maps made without 
any corrections. The amount of the shift was determined 
by aligning the optically thin part of the jet based on a 
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Figure 4. Plot of the positions for the inner jet components. Different symbols denote different frequencies with the filled symbols 
representing the mean positions at 22 (filled circles), 43 (filled squares) and 86 GHz (filled triangles). To prevent confusion, only the 
averaged mean positions are shown at 86 GHz. A clear frequency-dependent displacement of the jet component positions can be seen 
and is indicated by red arrows. Spectra are shown for those components, which are visible at more than 2 frequencies (insets). 



Table 3. Position shift of jet components. Listed are the component Id., the time-averaged value of 10 d of components' core separation 
(r) and PA (6) at 22, 43 and 86 GHz. The last four columns summarize the magnitude and direction (in degrees east of north) of the 
position shift between the two frequency pairs (22-43 GHz and 43—86 GHz) for each component. 



Id. 



22 GHz 
r 6 

(mas) (°) 



43 GHz 
r 6 

(mas) (°) 



86 GHz 

r e 

(mas) (° 



Ar 2 2-43 
(mas) 



W22-43 
(°) 



Ar 4 3_ s6 
(mas) 



043- 86 
(°) 



n ... ... 0.38±0.01 -3.3 ± 1.4 0.37±0.01 -6.8 ±1.1 

d 1.17±0.02 -10.0 ±0.9 1.34±0.01 -9.8 ± 0.4 1.43±0.03 -7.2 ± 1.4 0.17 ±0.02 -8.5 ± 7.1 

e 1.48±0.02 11.0±0.8 1.60±0.02 9.9±0.7 ... ... 0.12 ± 0.03 -3.7 ±13.5 

f 3.27±0.02 27.8±0.4 3.41±0.03 25.1±0.5 ... ... 0.21 ± 0.04 -21.7 ±10.0 



0.03±0.01 
0.11 ±0.03 



-118.6 ±30.0 
26.3 ±18.3 



two-dimensional cross-correlation technique, a method very 
similar to that described by Croke & Gabuzda (2008). In 
our case, we restored both images with a circular beam of 
0.5 mas. The derived shift of the 43 GHz-image relative to 
the 22 GHz-image was 0.02 mas to the east and —0.16 mas 
to the north (Ar22-43 = 0.16 mas, 622-43 = —7.1°), in a di- 
rection very close to the inner jet (PA ~ —10°). These shifts 
are consistent with the results obtained from model fitting 
(average shift of Ar22 - 43 = 0.16 mas, 622-43 = —11.5° for 
components d, e and f), indicating an overall relative shift 



effect between the two images. The derived shifts for other 
selected epochs also agree well with these values. 

A well known effect which could cause the frequency de- 
pendence of jet component position is the opacity-induced 
position shift of the VLBI core ('core shift', Marcaide & 
Shapiro 1984; Lobanov 1998) in a synchrotron-self absorbed 
jet (Blandford & Konigl 1979). The VLBI core position, 
which is defined by the r = 1 surface at the jet base, de- 
pends on the frequency, with the distance to the central 
engine, r corc cc u^ 1 ^ 1 , where fc r is related to the electron en- 
ergy distribution, the magnetic field, and the particle num- 



High-frequency very long baseline interferometry studies of NRAO 530 7 




J i i i i i i i i i L 

I 



Core separation (mas) 

Figure 5. Slice along PA = —10° through the total intensity 
maps of the inner 2 mas of NRAO 530 at 22 and 43 GHz on 2007 
May 15 and at 86 GHz on 2007 May 16. An increase of the core 
separation of components n and d with frequency can be inter- 
preted as opacity shift of the core position (see text). 

ber density (Konigl 1981). For a self-absorbed core which is 
in energy equipartition, one can show that k r — 1 (Lobanov 
1998). Recent measurements indeed show a good agreement 
with fc r = 1 (e.g. O'Sullivan & Gabuzda 2009). Between two 
frequencies [yi > vi), the core shift can be expressed as 

Ar C ore = (rcore— fcore) OC (v7 1/kl -vZ 1/kl ) |fe r= i(X — -.(3) 

1*2 1/1 

Since the absolute position information for VLBI observa- 
tions is lost due to the phase self-calibration performed in 
the imaging, core shifts are normally measured by referenc- 
ing the core to the optically thin jet components, whose po- 
sitions are expected to be frequency-independent. We used 
the optically thin component d (which is visible at all fre- 
quencies between 15 and 86 GHz, see Table 9) to measure 
the VLBI core positions shift between 86 GHz and the lower 
frequencies. We estimated k r = 0.78±0.29, consistent within 
error with k r = 1 as expected if the region in question is in 
equipartition (Fig. 7). 

The number of extragalactic radio sources with reported 
core shift has been increasing (Cai et al. 2007; Kovalev et al. 
2008; O'Sullivan & Gabuzda 2009, and references therein). 
Kovalev et al. (2008) studied a core shift sample of 29 AGN 
between 2.3 and 8.4 GHz. They found the magnitude of mea- 
sured core shifts between 2.3 and 8.4 GHz to range from —0.1 
to 1.4 mas with a median value of 0.44 mas 2 . For higher fre- 
quency pairs, we would expect that the core shift decreases 
since opacity effects become less important (equation 3). For 
a few sources where core shifts between high frequency pairs 
are available (O'Sullivan & Gabuzda 2009), the reported 
core shifts between 22 and 43 GHz (0.02-0.04 mas) are sig- 
nificantly smaller than those measured here for components 
d, e, and f (0.12-0.17 mas). In addition, the displacement 
of the jet components' position are not all in the same di- 
rection. In other words, we observe non-radial shifts, which 
may indicate two dimensional gradients of the jet parameters 

2 Positive values mean that the separation of the jet component 
measured at higher frequency is larger than the one measured at 
lower frequency in one frequency pair. 
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Figure 6. Spectral index map of NRAO 530 with a single contour 
(5 m Jy beam" 1 ) of the 22 GHz Stokes I image superimposed. The 
map was made with 22 and 43 GHz data on 2007 May 15 after 
applying for the 43 GHz map a relative shift of 0.02 mas to the 
east and -0.16 mas to the north (see text for details). Note that 
the oblique boundary between the regions of optically thick and 
thin emission near the core is just perpendicular to the inner jet 
direction along a PA = —10°. 




_ , 1 

20 40 60 80 100 

Frequency (GHz) 

Figure 7. Position shift plot for component d with 86 GHz as 
reference frequency. 

(magnetic field, particle number density), an effect which is 
expected to occur if a bent jet is observed at a small viewing 
angle. 

3.3 Flux density and structure variability on daily 
time-scales 

Our data also permit the study of flux density and struc- 
ture variability in the VLBI jet of NRAO 530 on a daily 
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time-scale. We show in Fig. 8 the variations of the flux den- 
sity of the model fit components with their modulation in- 
dices summarized in Table 4. On daily time-scales, the flux 
density of most model fit components shows no significant 
variations. Component g seems to be more variable (m = 53 
per cent). However, it is the weakest component represent- 
ing a diffuse region with flux density of only ~ 30 mjy This 
may in turn reflect our modelling uncertainties. Through a 
X 2 -test (see Table 4) we found the probability for variability 
for most of the components to be low. For component c and 
component n, the x 2 -test indicates the presence of interday 
variability with probability of 85 per cent at 22 GHz and 92 
per cent at 43 GHz for component c and 99.9 per cent at 
43 GHz for component n. We therefore conclude that the in- 
ner jet components are more variable than those further out. 
However, since component n is located very close to compo- 
nent c, the apparent variability may be caused by variations 
in how the flux of component c and n is divided up due 
for example to different uv coverage. An anti-correlation be- 
tween the fluxes of component c and n at 43 GHz (correlation 
coefficient: —0.75) seems to support this scenario, while it is 
less clear at 86 GHz (correlation coefficient: — 0.2f). Future 
VLBI studies are needed to confirm this effect. 

Similarly, the variations of the core separation and PA 
of the jet components are displayed in Fig. 9. The mod- 
ulation indices m, xt, and the probability for the studied 
source parameters to be variable are tabulated in Tables 5 
and 6. For the core separation, the scatter in component g 
is again much larger than what we expected from statistical 
errors. In addition, the large scatter for the core separation 
of component n and d at 86 GHz also indicates the presence 
of additional errors that are not accounted for by the general 
accepted 'rule of thumb', i.e. one-fifth beam size. This may 
be due to the sparse uv coverage at 86 GHz and the weakness 
of the components. For the PA, all the components show 
stationarity except component d at 86 GHz, which shows 
similar large scatter to the core separation. The limited uv 
coverage and calibration accuracy at 86 GHz, however, may 
cause artificial variability, making it necessary to confirm 
this effect in future VLBI studies. 

3.4 Jet kinematics at 15 GHz 

In this section, we investigate the kinematics of the jet com- 
ponents at 15 GHz. The plot of the core separation of the 
jet components as a function of time, as shown in Fig. 10, 
served as a tool for the cross-identification across epochs. 
We, however, note that the identification of the innermost 
jet component d at epochs before 2001 may be subject to 
ambiguity due to the gaps during 2000 and 2002. The sys- 
tematic change of flux density and PA (see later in this sec- 
tion), however, supports our identification scheme also for 
component d. Linear regression fits to the data after 2002 
yielded for the epoch of ejection a date, coinciding with the 
time of a radio flare (~ 2002.6, Hong et al. 2008; Aller et 
al. 2009). This may indicate that component d was newly 
ejected during that flare. Component e appears to separate 
from component f after 2003. This scenario is supported by 
a sudden change in the flux density of component f (whose 
flux density would be maintained after the separation by 
adding the flux density of component e) , the similar PAs of 
f and e in 2004, and the lack of components between d and 



fin earlier epochs before 2004. The most reliable identifica- 
tion corresponds to outer components f, g, h, i, j, and the 
newly ejected component n'. 

In order to derive estimates for the separation speeds of 
the individual components, we applied linear regression fits 
[f(t) = a ■ (t — to), where to is the component ejection time] 
to the separation of the components from the core versus 
the epoch of observations. The fitting results are shown in 
Fig. 10 as dashed lines. The angular separation rates, the 
time of ejection, and the apparent speeds are summarized 
in Table 7. 

Component n seems to show the maximum apparent 
speed (/3 apP = 26.5 ± 12.5). However, the speed had to be 
calculated on the basis of only two epochs and is therefore 
subject to large uncertainties. We note that component n' 
apparently was ejected around 2008.6, which coincides with 
the second moderate flare seen in late 2008 at millimeter 
wavelengths 3 and in early 2009 at centimeter wavelengths 
(Aller et al. 2009). Component f and g can be identified with 
components B and E seen by Feng et al. (2006). The derived 
apparent speeds (14 c and 21 c) are slightly faster than their 
findings (10c and 14c), based on their shorter time cover- 
age (~ 4 yr). The ejection time of both components lies in 
the range of the dramatic outburst phase between 1994 and 
1998. However, due to the expansion effects, the uncertainty 
of core separation for component g gets systematically larger 
at later epochs and, at the same time, it shows an accelera- 
tion indicating that a linear fit to its motion is only a rough 
approximation. Component e and the three outer compo- 
nents (h, i and j) reveal relatively slow motions with P^pp in 
the range of 2-9 c. 

The measured apparent velocities can be used to esti- 
mate the parameters of the jet components according to the 
superluminal motion equation 



where /3 apP and /3 are the apparent and the true velocity 
in units of speed of light c and 9 is the angle between the 
direction of motion and the line of sight. Making use of the 
apparent speed of component f (/3 ap p = 14.1 c), which is the 
most reliably determined fastest component, w e can set a 
lower limit on the Lorentz factor V, r m i n = yl + $*pp = 
14.1, and derive maximum allowable viewing angle (? max of 
8?1. The critical viewing angle, which maximizes the /? a pp 
for a given /3, is # C ri = arcsin(i). For component f, 9 CTl is 
4?1. Correspondingly, we can also derive a Doppler factor 
5mtn using the minimum Lorentz factor F m i n and # cr j. For 
component f, 5 m in is 14.1. For even smaller angles (6 —¥ 0), 
8 tends to the limit of ~ 2r = 28.2. 

Fig. 11 shows the relation between the apparent speed, 
the Doppler factor and the viewing angle. For a given 
Lorentz factor, the apparent speed first shows a shallow 
increase with decreasing viewing angle and then decreases 
sharply towards very small angles. Therefore, equation (4) 
has, for a given apparent speed, two solutions for the view- 
ing angle. However, the Doppler factor, which has a strong 
influence on the observed flux density due to the beaming ef- 
fects, shows a different orientation dependence and increases 

3 cf. http:/ /smal.sma.hawaii.edu/callist/callist.html?data=1733- 
130 
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Figure 8. Variations of the flux density of model fit components on daily time-scales. Left: outer jet components (i, g, f and e), Right: 
inner components (d, n and c). Note that components i and g are only visible at 22 GHz. 



Table 4. Intcrday flux variability characteristics of model fit components of NRAO 530. The columns list the component Id., the 
modulation index, reduced xi an d probability for the observable being variable. 
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Figure 9. Variations of the core separation (left) and PA (right) of model fit components on daily time-scales. Different symbols denote 
different frequencies (circles: 22 GHz; squares: 43 GHz; triangles: 86 GHz). The core separation of component i and e has been shifted for 
display purposes. 
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Table 5. Interday variability characteristics of core separation for the jet components. The columns list the component Id., the modulation 
index, reduced xl an d probability for the observable being variable. 
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Table 6. Interday PA variability characteristics of model fit components for NRAO530. The columns list the component Id., the 
modulation index, reduced xl an d probability for the observable being variable. 
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Figure 10. Separation from the core as a function of time for the VLBI jet components seen in NRAO 530. The core separations of 
component j have been shifted 13 mas downwards for display purposes, and the component x is not shown. The dashed lines are linear 
regression fits to the data. 
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Table 7. Results of the linear fits for the component motion at 
15 GHz. Listed are the component Id., number of epochs for the 
fit, proper motion fi, apparent velocity /3 app and ejection time of 
each component tg. 



W. # jU /3app t 

(mas yr" 1 ] (c) (yr) 



n' 


2 


0.55±0.26 


26.5±12.5 


2008.6±0.4 


cl 


14 


0.26±0.02 


12.4±1.0 


2002.5±0.3 


c 


11 


0.13±0.02 


6.2±1.1 


1994.4±2.4 


f 


17 


0.29±0.01 


14.1±0.5 


1996.0±0.4 


g 


11 


0.44±0.03 


21.3±1.5 


1996.4±0.5 


h 


9 


0.05±0.05 


2.3±2.5 


1876.6±133.6 a 


i 


14 


0.19±0.01 


9.2±0.5 


1960.8±3.6 


j 


17 


0.13±0.01 


6.3±0.5 


1823.1±14.7 



a From formal back-extrapolation 
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Figure 11. Apparent speed (upper panel) and Doppler factor 
(lower panel) as a function of viewing angle for component f with 
a Lorentz factor of 14.1. The vertical line indicates the critical 
angle (4? 1) at which the apparent speed is maximized. 

continuously with decreasing viewing angle. Based on these 
signatures and how the apparent speed and flux density of 
the jet component evolves, we can estimate the jet geome- 
try. For component g, we found that the apparent jet speed 
increases and the flux density decreases at later epochs [see 
Figs 10 and 12 (right)], indicating that the viewing angle 
takes the small viewing angle solution and is approaching the 
critical value. These signatures are reminiscent of a three- 
dimensional (spatially) bent jet. 

In Fig. 12 (left), we present the evolution of the PA of 
jet components with time. One can clearly see a systematic 
variation of the PAs for components d, e, h and i is in con- 
trast to the stationarity of the components f, g and j. The 
newly ejected component n also shows a similar position 
angle rotation (not shown). The changes of PA with time 



clearly indicate a non-radial, non-ballistic motion of the jet 
components. 

In Fig. 12 (right), the flux density evolution of the mod- 
elled components at 15 GHz is displayed. Interestingly, the 
flux density evolution of the modelled components is charac- 
terized by several distinct behaviors. For components i and 
j, the flux density remains nearly constant, while for compo- 
nents g and h, the flux density continuously decreases with 
time. It seems that the drop in the light curve of component 
f around 2004 was due to the separation of e from it. Com- 
ponent e, however, shows a clearly increasing trend. Com- 
ponent d shows convolved behavior before the peak around 
2004 and fades afterwards as it travels down the jet. The 
core component c shows two intermediate flares, which co- 
incide with the ejection of component xi and probably also 
with the ejection of component d, if it was ejected around 
2002 according to the straightforward back-extrapolation. 
Component c and f show correlated variations before 2005 
(significant at the 90 per cent confidence level), and the 
flux changes in components c and e are also correlated after 
~ 2004, significant at the level of 95 per cent, for which the 
physical reason is unclear, but could be due, for example, to 
differential Doppler boosting. 



3.5 Morphology and its evolution 

A comparison of the projected trajectories of jet components 
in 2007 and 1997 (Feng et al. 2006) indicated that the mor- 
phology of the jet changed significantly. It is obvious that 
the major differences are in the regions close to the core, 
in contrast to the relative 'position' stationarity of the out- 
ermost component. This can be attributed to the relatively 
fast motion of the inner jet components. Fig. 13 (left) shows 
the evolution of the jet ridge line over the last 10 years at 
15 GHz. This plot illustrates that the morphology of the jet 
changes from epoch to epoch due to the motion of the jet 
components. Shown in Fig. 13 (right) is the evolution of the 
inner jet axis using component d and e as tracers. A jet PA 
swing of ~ 3? 4 per year can be estimated for both compo- 
nents between 2004 and 2009. For component d, the average 
swing speed is also ~ 3?4 per year between 1999 and 2009. 
This indicates that the entire inner jet is swing together. 



4 DISCUSSION 

In this section, we discuss the implications of the jet bend- 
ing and swing which we observed for the pc-scale jet in 
NRAO 530. Recently, an increasing number of blazars show 
clear evidence of jet PA swing, either periodic or erratic, e.g. 
3C 120 (Gomez et al. 1999), BL Lac (Stirling et al. 2003), 
OJ 287 (Tateyama & Kingham 2004), 3C 273 (Savolainen et 
al. 2006) and NRAO 150 (Agudo et al. 2007). The change of 
the direction of the inner jet axis may be also responsible for 
the jet curvature observed at pc and kpc scales (e.g. Gomez 
et al. 1999). Since the jet axis is close to the observer's line 
of sight (as indicated by the observed apparent superluminal 
motion), it is certain that the apparently curved morphology 
results from amplification effects by projection. However, the 
physical reason for the intrinsic jet bending and swing is still 
unclear. 
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Figure 12. Time evolution of PA for components d-j (left) and 15 GHz flux density of each model fit component (right). 
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Figure 13. Evolution of the projected jet axis. Left: the jet ridge line between 1999 and 2009. Right: motion of inner jet component d 
and e in rectangular coordinates between 2004 and 2009. Close epochs arc not shown in order to avoid cluttering the plot. 



Changes in the PA of the jet can be attributed to pre- 
cession effects. Jet precession could arise from scenarios re- 
lying on a single rotating black hole (e.g. Caproni, Mosquera 
Cuesta & Abraham 2004) or a binary black hole system (e.g. 
Begelman, Blandfordv & Rees 1980; Lobanov & Roland 
2005; Roland et al. 2008; Britzen et al. 2010). (Quasi)- 
periodicities in the radio and optical light curves and helical 
motions of jet components argue in favor of precession mod- 
els (e.g. Kudryavtseva et al. 2011; Liu et al. 2010). Further- 
more, regular variations of the component ejection PA with 
a likely period of a few years were observed in a number of 
sources, e.g. 3C 273 (Abraham & Romero 1999). Abraham & 
Romero (1999) proposed a simple ballistic precessing model 
of the jet 'nozzle' to explain the kinematics in the quasar 
3C 273. Similar models have been proposed and applied to 
other objects, e.g. 3C 279 (Abraham & Carrara 1998), 3C 
345 (Caproni & Abraham 2004), and OJ 287 (Tateyama & 
Kingham 2004). In the case of NRAO 530, however, the non- 



ballistic motion of the components suggests that a ballistic 
treatment of the jet is probably not appropriate. 

Another alternative explanation are plasma instabilities 
in the context of general relativistic magnetohydrodynami- 
cal (GRMHD) jet formation and acceleration (Blandford & 
Znajek 1977; Blandford & Payne 1982), most importantly, 
Kelvin-Helmholtz (KH) instability (e.g. Hardee & Norman 
1988; Perucho et al. 2006). KH instabilities can be driven 
by velocity shear at the boundary layers between the jet 
and the ambient medium. Ideas of current-driven (CD) kink 
instability are also considered (e.g. Mizuno et al. 2009). It 
is shown that the strength of CD kink instability strongly 
depends on the magnetic pitch profile. However, it is still 
not clear whether velocity shear or current driven instabil- 
ity are responsible for the observed bending structure and 
jet swing, or whether they act at different spatial scales. 
Furthermore, instabilities may only determine the morphol- 
ogy and kinematics for jet regions at large scales (beyond 
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~100 pc), while the inner part of the flow is strongly shock- 
dominated (Lobanov & Zensus 1999; Lobanov 2010). For 
the jet features in NRAO 530, the non-ballistic motion is 
observed in the inner few pc, a region where the growth of 
instabilities may be weakened by the magnetic field, similar 
to the case of NRAO 150 (Agudo et al. 2007). 



5 CONCLUSIONS 

We studied the pc-scale jet properties of NRAO 530 by us- 
ing high-resolution multi-epoch VLBI observations at 15, 
22, 43 and 86 GHz. We identified component c as the core 
based on its unambiguous spectral turnover, compactness 
and variability using the quasi-simultaneous VLBI data ob- 
tained in 2007. The frequency-dependent position shifts are 
found to be significant for three components (component d, 
e and f) with mean radial shift of 0.15 ± 0.02 mas between 
22 and 43 GHz. The position shift can be interpreted as be- 
ing due to a core shift effect in a synchrotron self-absorbed 
jet. Interestingly, we find, for the first time, that jet compo- 
nents show a two-dimensional position shift. This indicates 
that the variations of the optical depth, which induces the 
core-shift, have also a non-radial dependence. Our 10-d mon- 
itoring observations allow us to investigate the jet variability 
characteristics on daily time-scales for NRAO 530. We find 
the inner jet to more variable than the outer jet. 

We studied the jet kinematics with individual compo- 
nents using the archival data at 15 GHz. We observed super- 
luminal motion for eight components with apparent speeds 
/9 apP in the range of 2-26c. From back-extrapolation, the 
ejection of component n' can be associated with the mod- 
erate flare in the radio light curve around 2009. However, 
the identification of the ejection time for another compo- 
nent d is subject to uncertainties, and so is the link between 
its ejection and the moderate flare around 2002. From the 
measured apparent jet speed, we set limits to the jet pa- 
rameters. For the reliably measured fastest component, we 
obtained a minimum Lorentz factor of 14 and Doppler fac- 
tors in the range of 14-28. Different jet components show 
different kinematical behavior and flux density evolution, 
which can be understood as a consequence of a rotating jet. 
Using the variations of the apparent speed and flux den- 
sity, we can estimate the jet geometry. The morphology of 
NRAO 530 was found to have changed significantly over the 
past 10 years between 1997 and 2007. The variations are 
more pronounced near the core in contrast to the appar- 
ent 'stationarity' of components further out. Moreover, we 
found significant variations of the jet PA and no common 
path existed for all the jet components. For components d 
and e , we estimated a mean swing of ~3?4 per year, which is 
partly responsible for the variations of the inner jet structure 
and indicates the swing of the entire inner jet. We discuss 
our results in the light of jet precession and models of jet 
instabilities. 

To summarize, we conclude that the jet of NRAO 530 
is far more complex than previously known. The observed 
non-ballistic motion and a pronounced change of the jet ori- 
entation with a rate of 3? 4 per year suggest that NRAO 530 
is another example of a "swinging" jet, similar to e.g. BL Lac 
(Stirling et al. 2003) and NRAO 150 (Agudo et al. 2007). 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the on- 
line version of this article: 

Fig. 14. Clean maps of NRAO 530 at 22, 43 and 86 GHz 
obtained during our 10-d observations in 2007 May. 

Fig. 15. Clean maps of NRAO 530 at 15 GHz from the 
archival data between 1999-2009. 

Table 8. Description of VLBA images of NRAO 530 
shown in Figs 14 and 15. 

Table 9. Model-fitting results including the epoch of ob- 
servation, the identification of components and the fitted 
parameters and their uncertainty. 
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Figure 14. Clean maps of NRAO 530 at 22, 43, and 86 GHz in May 2007. 
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Figure 14. -continued. 
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Figure 14. -continued. 
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Table 8: Description of VLBA images of NRAO 530 shown in Fig. 14 and Fig. 15. The columns list the epoch of observation, 
the observing frequency, the peak flux density, the parameters of the restoring elliptical Gaussian beam: the full width at 
half maximum (FWHM) of the major and minor axes and the position angle of the major axis, and the contour levels of the 
image, expressed as a percentage of the peak intensity. 



Restoring Beam 



Epoch 


V 


•Speak 


Major 


Minor 


P.A. 


Contours 




GHz 


Jy/beam 


mas 


mas 


deg 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


2007.370 


22 


1.32 


1.19 


0.344 


-12.8 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.33 


0.537 


0.193 


-7.7 


-0.3, 0.3, 0.6, ... 


76.8 




86 


0.74 


0.2 


0.2 





-1, 1, 2, 


.., 64 


2007.373 


22 


1.31 


1.29 


0.344 


-15.0 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.33 


0.628 


0.19 


-12.9 


-0.3, 0.3, 0.6, ... 


76.8 




86 


1.20 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.375 


22 


1.29 


1.14 


0.367 


-9.2 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.41 


0.794 


0.252 


-11.2 


-0.3, 0.3, 0.6, ... 


76.8 




86 


1.06 


0.2 


0.2 





-1, 1, 2, 


.., 64 


2007.378 


22 


1.38 


1.24 


0.367 


-10.0 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.44 


0.72 


0.314 


-12.0 


-0.3, 0.3, 0.6, ... 


76.8 




86 


1.25 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.381 


22 


1.32 


1.15 


0.374 


-10.1 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.44 


0.83 


0.319 


-10.0 


-0.3, 0.3, 0.6, ... 


76.8 




86 


1.04 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.384 


22 


1.27 


1.12 


0.42 


-6.88 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.37 


0.571 


0.186 


-8.6 


0.3, 0.6, 1.2, ... 


76.8 




86 


1.13 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.386 


22 


1.31 


1.09 


0.403 


-5.95 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.37 


0.604 


0.198 


-8.03 


0.3, 0.6, 1.2, ... 


76.8 




86 


1.01 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.389 


22 


1.37 


1.26 


0.424 


-8.68 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.41 


0.628 


0.199 


-9.78 


0.3, 0.6, 1.2, ... 


76.8 




86 


1.04 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.392 


22 


1.34 


1.20 


0.43 


-6.14 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.37 


0.553 


0.197 


-7.02 


0.3, 0.6, 1.2, ... 


76.8 




86 


1.19 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


2007.395 


22 


1.33 


1.37 


0.405 


-11.8 


-0.3, 0.3, 0.6, ... 


76.8 




43 


1.45 


0.663 


0.183 


-12.8 


0.3, 0.6, 1.2, ... 


76.8 




86 


1.13 


0.2 


0.2 





-0.5, 0.5, 1, 


.., 64 


1999.712 


15 


2.07 


1.35 


0.495 


-6.8 


0.25, 0.5, 1, 


.., 64 


2000.030 


15 


1.83 


1.42 


0.505 


-9.0 


-0.25, 0.25 ,0.5, 


.., 64 


2000.333 


15 


1.99 


1.35 


0.509 


-2.6 


-0.3, 0.3 ,0.6, ... 


76.8 


2002.773 


15 


4.73 


1.51 


0.547 


-7.7 


0.1, 0.2 ,0.4, ... 


51.2 


2003.258 


15 


4.15 


1.34 


0.55 


-2.2 


0.15, 0.15 ,0.3, ... 


76.8 


2003.337 


15 


4.23 


1.36 


0.544 


-2.8 


0.15, 0.3 ,0.6, ... 


76.8 


2004.115 


15 


3.36 


1.35 


0.519 


-5.2 


0.15, 0.3 ,0.6, ... 


76.8 


2005.225 


15 


2.81 


1.2 


0.483 


-3.9 


-0.2, 0.2 ,0.4, ... 


51.2 


2005.562 


15 


2.29 


1.37 


0.529 


-6.0 


0.2, 0.4 ,0.8, ... 


51.2 


2006.515 


15 


1.72 


1.38 


0.524 


-8.3 


-0.25, 0.25, 0.5, 


.., 64 


2007.099 


15 


1.54 


1.59 


0.551 


-9.3 


-0.3, 0.3, 0.6, ... 


76.8 


2007.441 


15 


1.20 


1.34 


0.568 


-5.7 


-0.35, 0.35, 0.7, ... 


89.6 


2008.333 


15 


2.76 


1.29 


0.509 


-4.2 


-0.15, 0.15, 0.3, ... 


76.8 


2008.544 


15 


2.91 


1.33 


0.524 


-3.7 


-0.15, 0.15, 0.3, ... 


76.8 


2008.757 


15 


3.19 


1.34 


0.559 


-4.2 


-0.2, 0.2 ,0.4, ... 


51.2 


2009.153 


15 


3.63 


1.31 


0.524 


-5.4 


0.15, 0.3, 0.6, ... 


76.8 


2009.482 


15 


3.19 


1.30 


0.51 


-4.5 


0.25, 0.5, 1, 


.., 64 
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Table 9: Model-fitting results for NRAO530. 



Epoch 


Id. 


Flux 


Core Separation 


P.A. 


Size 






[mJy] 


[mas] 


[degree] 


[mas] 




(I) V 


= 22 GHz 




2007.370 


c 


1332.9±23.6 


O.OOiO.OO 


O.OiO.O 


0.08±0.01 




d 


385.9±40.5 


1.10±0.06 


-9.5±3.1 


0.28±0.02 




c 


377.8±72.1 


1.52±0.08 


9.6±2.9 


0.87±0.15 




f 


99.1±26.5 


3.41±0.10 


26.1±1.6 


0.81±0.20 




g 


29.5±11.2 


4.52±0.22 


7.0±2.8 


0.32±0.09 




i 


138.3±44.2 


8.75±0.30 


13.4±2.0 


1.99±0.61 


2007.373 


c 


1317.6±60.3 


O.OOiO.OO 


O.OiO.O 


0.10±0.01 




d 


390.3±45.4 


1.16±0.06 


-9.6±3.0 


0.28±0.03 




c 


358.2±74.2 


1.48±0.08 


10.8±3.0 


0.83±0.16 




f 


76.1±19.8 


3.28±0.06 


28.1±1.1 


0.56±0.12 




g 


52.6±23.6 


4.56±0.31 


9.4±3.8 


1.49±0.61 




i 


116.7±38.6 


8.94±0.26 


13.6±1.6 


1.64±0.52 


2007.375 


c 


1289.7±43.2 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




d 


435.9±54.6 


1.15±0.06 


-9.6±3.0 


0.32±0.03 




c 


298.9±58.2 


1.54±0.07 


13.2±2.5 


0.78±0.14 




f 


98.8±24.5 


3.27±0.10 


26.9±1.8 


0.88±0.20 




g 


26.3±16.2 


4.82±0.20 


6.4±2.4 


0.41±0.19 




i 


123.0±42.1 


8.88±0.28 


13.0±1.8 


1.74±0.57 


2007.378 


c 


1387.8±30.8 


O.OOiO.OO 


O.OiO.O 


0.14±0.01 




(1 


413.1±36.2 


1.17±0.06 


-10.7±2.9 


0.25±0.02 




c 


345.6±63.1 


1.54±0.06 


10.6±2.2 


0.74±0.12 




f 


97.8±20.9 


3.13±0.06 


29.6±1.1 


0.48±0.08 




g 


14.2±13.2 


5.19±0.16 


7.3±1.7 


0.54±0.31 




i 


83.3±29.1 


8.87±0.21 


15.1±1.3 


1.29±0.42 


2007.381 


c 


1313.5±27.7 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




(1 


392.4±25.0 


1.18±0.06 


-9.9±2.9 


0.27±0.01 




c 


352.3±62.8 


1.42±0.07 


10.8±2.7 


0.82±0.13 




f 


89.0±16.9 


3.31±0.06 


27.3±1.0 


0.70±0.12 




g 


58.3±27.6 


4.31±0.39 


8.4±5.2 


1.78±0.79 




i 


121.1±35.8 


8.89±0.20 


13.4±1.6 


1.73±0.49 


2007.384 


c 


1249.7±28.4 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




d 


388.1±22.8 


1.16±0.06 


-10.7±3.0 


0.30±0.01 




c 


364.1±52.5 


1.43±0.06 


11.2±2.4 


0.85±0.10 




f 


94.1±23.0 


3.33±0.08 


26.6±1.3 


0.73±0.15 




g 


22.0±16.7 


4.72±0.20 


6.5±2.4 


0.35±0.20 




i 


100.5±30.1 


8.95±0.20 


14.3±1.3 


1.43±0.39 


2007.386 


c 


1298.3±42.6 


O.OOiO.OO 


O.OiO.O 


0.08±0.01 




d 


409.7±47.1 


1.19±0.06 


-9.4±2.9 


0.32±0.03 




c 


302.6±57.5 


1.49±0.07 


11.6±2.7 


0.83±0.14 




f 


89.7±23.6 


3.35±0.10 


26.6±1.7 


0.82±0.19 




g 


34.3±12.7 


4.58±0.26 


7.0±3.2 


0.42±0.12 




i 


112.0±36.9 


8.91±0.26 


13.6±1.7 


1.66±0.52 


2007.389 


c 


1336.8±41.5 


O.OOiO.OO 


O.OiO.O 


0.06±0.01 




d 


403.2±20.2 


1.22±0.06 


-9.8±2.8 


0.24±0.01 




c 


338.1±54.8 


1.45±0.06 


11.8±2.4 


0.74±0.11 




f 


82.6±11.3 


3.22±0.06 


27.9±1.1 


0.55±0.06 




g 


80.3±30.5 


4.03±0.36 


10.9±5.1 


2.02±0.72 




i 


118.5±29.1 


9.04±0.18 


13.0±1.1 


1.56±0.36 


2007.392 


c 


1344.6±34.4 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




d 


399.8±31.5 


1.22±0.06 


-10.4±2.8 


0.26±0.02 




c 


328.5±40.2 


1.48±0.06 


11.4±2.3 


0.74±0.08 




f 


79.9±15.2 


3.24±0.06 


27.9±1.1 


0.56±0.09 




g 


79.1±24.5 


3.86±0.24 


10.5±3.6 


1.73±0.49 




i 


140.2±30.5 


8.88±0.19 


12.9±1.2 


1.83±0.37 


2007.395 


c 


1294.6±37.0 


O.OOiO.OO 


O.OiO.O 


0.08±0.01 




d 


364.7±41.3 


1.15±0.06 


-10.0±3.0 


0.25±0.02 
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Table 9 — continued 


Epoch 


Id. 


Flux 


Core Separation 


P.A. 


Size 






[mJy] 


[mas] 


[degree] 


[mas] 




e 


416.2±75.9 


1.48±0.08 


7.7±3.2 


0.98±0.17 




f 


89.5±21.9 


3.35±0.07 


28.7±1.1 


0.63±0.13 




g 


39.9±19.2 


4.37±0.17 


8.1±2.2 


0.83±0.34 




i 


118.9±34.3 


8.92±0.22 


12.8±0.4 


1.63±0.45 




(I) V 


= 43 GHz 




2007.370 


c 


1258.0±60.7 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




n 


272.9±22.4 


0.38±0.03 


-3.2±4.5 


0.08±0.01 




(1 


221.2±27.7 


1.31±0.03 


-10.4±1.3 


0.21±0.02 




e 


250.3±40.2 


1.56±0.05 


7.5±2.0 


0.70±0.11 




f 


56.1±20.2 


3.44±0.12 


24.4±2.0 


0.68±0.24 


2007.373 


c 


1205.8±52.9 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




n 


295.5±17.4 


0.35±0.02 


-3.4±4.9 


0.10±0.01 




(1 


217.7±42.0 


1.32±0.03 


-10.2±1.3 


0.22±0.04 




c 


245.8±49.2 


1.57±0.07 


8.6±2.5 


0.72±0.14 




f 


48.9±16.8 


3.50±0.11 


24.9±1.8 


0.70±0.23 


2007.375 


c 


1303.1±47.9 


O.OOiO.OO 


O.OiO.O 


0.01±0.01 




11 


233.1±19.2 


0.42±0.03 


-5.1±4.1 


0.12±0.01 




(1 


235.9±32.8 


1.34±0.03 


-9.4±1.3 


0.23±0.02 




e 


195.5±35.4 


1.67±0.06 


11.1±1.9 


0.67±0.11 




f 


56.8±14.1 


3.47±0.07 


25.6±1.1 


0.60±0.14 


2007.378 


c 


1303.1±29.8 


O.OOiO.OO 


O.OiO.O 


0.02±0.01 




n 


261.9±22.0 


0.39±0.03 


-2.9±4.4 


O.lliO.Ol 




(1 


276.4±40.4 


1.36±0.03 


-8.3±1.3 


0.27±0.03 




c 


181.0±41.3 


1.70±0.07 


12.5±2.3 


0.66±0.14 




f 


55.6±20.8 


3.45±0.13 


25.2±2.1 


0.73±0.25 


2007.381 


c 


1269.8±33.2 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




11 


302.8±19.1 


0.36±0.03 


-4.0±4.8 


0.12±0.01 




(1 


261.9±36.6 


1.36±0.03 


-8.5±1.3 


0.25±0.02 




c 


181.1±35.6 


1.67±0.06 


11.9±2.0 


0.66±0.11 




f 


57.0±13.2 


3.50±0.07 


25.0±1.2 


0.67±0.14 


2007.384 


c 


1274.3±18.2 


O.OOiO.OO 


O.OiO.O 


0.02±0.01 




11 


234.4±22.0 


0.37±0.03 


-3.0±4.6 


0.05±0.01 




(1 


206.9±25.7 


1.33±0.03 


-10.5±1.3 


0.20±0.02 




c 


253.5±57.9 


1.56±0.08 


8.0±2.9 


0.73±0.15 




f 


51.9±17.7 


3.33±0.10 


25.0±1.8 


0.65±0.20 


2007.386 


c 


1287.2±51.5 


O.OOiO.OO 


O.OiO.O 


0.01±0.01 




11 


251.1±16.2 


0.36±0.03 


-3.4±4.8 


0.06±0.01 




(1 


239.4±32.1 


1.34±0.03 


-10.0±1.3 


0.22±0.02 




c 


226.2±55.1 


1.64±0.08 


9.9±2.9 


0.71±0.16 




f 


52.3±15.7 


3.34±0.10 


25.5±1.7 


0.70±0.20 


2007.389 


c 


1265.3±11.6 


O.OOiO.OO 


O.OiO.O 


0.02±0.01 




11 


241.7±21.8 


0.37±0.03 


-3.2±4.6 


0.06±0.01 




(1 


220.1±30.2 


1.34±0.03 


-10.1±1.3 


0.21±0.02 




c 


246.8±58.8 


1.57±0.09 


9.0±3.1 


0.75±0.17 




f 


52.8±14.9 


3.18±0.08 


24.6±1.5 


0.63±0.17 


2007.392 


c 


1308.5±18.9 


O.OOiO.OO 


O.OiO.O 


0.02±0.01 




11 


214.4±29.7 


0.39±0.03 


-2.7±4.4 


0.02±0.01 




(1 


212.8±27.6 


1.33±0.03 


-10.2±1.3 


0.20±0.02 




c 


243.7±38.4 


1.61±0.06 


8.7±2.0 


0.76±0.11 




f 


53.3±17.0 


3.42±0.11 


25.6±1.9 


0.76±0.23 


2007.395 


c 


1337.4±19.0 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




11 


210.5±13.4 


0.37±0.03 


-2.3±4.6 


0.02±0.01 




(1 


239.1±24.6 


1.34±0.03 


-10.5±1.3 


0.20±0.02 




c 


245.5±32.9 


1.54±0.05 


10.2±2.0 


0.83±0.11 




f 


50.5±14.5 


3.43±0.09 


24.6±1.5 


0.67±0.18 




(I) V 


= 86 GHz 




2007.370 


c 


1214.4±186.3 


O.OOiO.OO 


0.0 ±0.0 


0.16±0.02 


2007.373 


c 


1187.5±106.8 


O.OOiO.OO 


0.0 ±0.0 


0.06±0.01 
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Table 9 — continued 



Epoch 


Id. 


Flux 


Core Separation 


P.A. 


Size 






[mJy] 


[mas] 


[degree] 


[mas] 




n 


73.1±37.0 


0.45±0.02 


-9.6±2.5 


0.08±0.01 




d 


107.4±21.4 


1.55±0.05 


-9.7±1.8 


0.08±0.02 


2007.375 


c 


1077.2±94.0 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




n 


87.2±39.3 


0.37±0.02 


-6.7±3.1 


0.06±0.01 




d 


116.6±49.0 


1.39±0.05 


-4.9±2.1 


0.08±0.02 


2007.378 


c 


1194.6±70.3 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




n 


96.9±37.8 


0.35±0.04 


-4.5±6.0 


0.12±0.07 




d 


109.8±56.2 


1.55±0.07 


-0.5±2.6 


0.30±0.14 


2007.381 


c 


1088.2±39.0 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




n 


59.4±11.4 


0.37±0.02 


-9.3±3.1 


0.07±0.01 




d 


114.6±52.4 


1.32±0.10 


-11.7±4.1 


0.44±0.19 


2007.384 


c 


1170.0±54.9 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




n 


72.2±39.2 


0.35±0.02 


-6.8i3.3 


0.03±0.01 




d 


53.9±27.1 


1.47±0.05 


-5.4±1.9 


0.06±0.02 


2007.386 


c 


1067.4±88.2 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




n 


71.0±36.8 


0.31±0.02 


-4.3i3.9 


0.09±0.04 




d 


53.4±24.5 


1.43±0.05 


-5.6±2.0 


0.15±0.05 


2007.389 


c 


1091.5±81.2 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




n 


125.9±49.7 


0.43±0.06 


-12.2±8.4 


0.35±0.13 




d 


86.3±30.4 


1.33±0.04 


-12.3±1.6 


0.24±0.08 


2007.392 


c 


1192.8±75.8 


O.OOiO.OO 


O.OiO.O 


O.OliO.Ol 




11 


55.2±29.4 


0.35±0.02 


-7.3i3.3 


0.04±0.01 




d 


142.5±40.7 


1.39±0.05 


-7.8±2.1 


0.21±0.05 


2007.395 


c 


1131.2±60.9 


O.OOiO.OO 


O.OiO.O 


O.OliO.Ol 




n 


75.7±29.3 


0.36±0.02 


-0.9i3.2 


0.03±0.01 




d 


116.1±47.9 


1.57±0.11 


4.8±3.9 


0.59±0.21 




(II) V 


= 15 GHz 




1999.712 


c 


2013.6±48.7 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




d 


297.6±13.9 


0.47±0.06 


36.8±7.3 


0.25±0.01 




f 


145.5±12.0 


1.10±0.06 


29.1±3.1 


0.09±0.01 




g 


483.2±55.8 


1.55±0.06 


6.3±2.2 


0.93±0.10 




X 


36.5±9.7 


3.00±0.06 


16. Oil. 1 


0.43i0.09 




h 


217.6±43.2 


5.82±0.20 


13.4±1.9 


2.03±0.39 




j 


108.9±46.3 


23.06±0.93 


1.0±2.3 


4.42±1.86 


2000.030 


c 


1679.3±31.3 


O.OOiO.OO 


O.OiO.O 


0.06±0.01 




d 


303.5±25.5 


0.37±0.06 


18.6±9.2 


0.45±0.03 




f 


228.7±27.8 


1.12±0.06 


31.0±3.1 


0.35±0.03 




g 


320.0±33.2 


1.88±0.06 


6.2±1.8 


0.78±0.07 




X 


22.9±6.8 


2.64±0.06 


26.4±1.3 


0.11±0.02 




h 


197.1±43.1 


5.95±0.24 


12.8±2.3 


2.30±0.48 




j 


98.2±43.5 


23.79±1.31 


0.6±3.1 


5.94±2.61 


2000.333 


c 


1691.1±51.2 


O.OOiO.OO 


O.OiO.O 


0.07±0.01 




d 


489.5±25.3 


0.32±0.06 


17.0±10.6 


0.38±0.02 




f 


275.1±27.9 


1.24±0.06 


29.7±2.8 


0.40±0.03 




g 


391.6±42.3 


2.01±0.06 


6.3±1.7 


0.85±0.08 




X 


19.3±8.3 


2.66±0.06 


24.5±1.3 


0.06±0.02 




h 


219.6±40.8 


5.95±0.19 


13.4±1.8 


2.12±0.38 




j 


119.6±43.5 


23.18±0.87 


0.6±2.2 


4.89±1.75 


2002.773 


c 


4150.2±81.5 


O.OOiO.OO 


O.OiO.O 


0.03±0.01 




d 


693.7±59.5 


0.25±0.06 


9.7±13.5 


0.20±0.01 




f 


390.5±41.1 


1.85±0.06 


29.4±1.9 


0.59±0.05 




g 


209.3±22.5 


2.45±0.06 


4.6±1.4 


0.85±0.08 




h 


129.8±23.0 


5.30±0.21 


8.7±2.3 


2.47±0.42 




i 


118.6±27.6 


7.77±0.15 


16.4±1.1 


1.44±0.31 




j 


128.1±33.3 


23.38±0.53 


0.2±1.3 


4.15±1.07 


2003.258 


c 


3870.0±88.9 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




d 


494.0±26.3 


0.34±0.06 


7.8±10.0 


0.48±0.02 




f 


374.0±36.6 


1.91±0.06 


29.8±1.8 


0.67±0.05 



High-frequency very long baseline inter ferometry studies of NRAO 530 



Table 9 — continued 



Epoch 


Id. 


Flux 


Core Separation 


P.A. 


Size 






[mjy] 


[mas] 


[degree] 


[mas] 




g 


198.5±31.8 


2.69±0.07 


5.9±1.5 


0.99±0.14 




h 


138.2±28.6 


6.06±0.25 


9.3±2.4 


2.51±0.50 




i 


118.2±25.6 


8.04±0.13 


17.0±0.9 


1.31±0.26 




j 


136.4±42.6 


23.34±0.58 


0.1±1.4 


3.78±1.16 


2003.337 


c 


3907.3±79.1 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




(1 


456.4±31.4 


0.36±0.06 


8.8±9.5 


0.21±0.01 




f 


348.7±38.3 


1.85±0.06 


30.8±1.9 


0.59±0.05 




g 


295.0±45.0 


2.53±0.09 


8.3±1.9 


1.23±0.17 




h 


146.0±25.9 


5.97±0.22 


9.4±2.2 


2.61±0.45 




i 


131.9±25.9 


8.00±0.11 


16.9±0.8 


1.28±0.23 




j 


176.1±51.8 


23.57±0.61 


0.2±1.5 


4.22±1.23 


2004.115 


c 


2190.6±39.0 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




(1 


1315.0±31.0 


0.27±0.06 


6.9±12.5 


O.lliO.Ol 




c 


145.1±17.2 


1.21±0.06 


25.3±2.8 


0.83±0.09 




f 


167.5±16.1 


2.42±0.06 


26.8±1.4 


0.66±0.05 




g 


107.9±19.9 


3.05±0.09 


3.6±1.6 


1.03±0.17 




ll 


96.5±19.1 


6.16±0.23 


9.1±2.1 


2.42i0.46 




i 


111.4±18.2 


8.28±0.09 


16.4±0.6 


1.22±0.18 




j 


108.1±31.0 


23.66±0.53 


0.3±1.3 


3.75±1.06 


2005.225 


c 


2290.5±74.5 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




(1 


874.2±48.5 


0.49±0.06 


2.0±7.0 


0.21±0.01 




c 


141.0±21.6 


1.50±0.06 


21.5±2.3 


0.69±0.09 




f 


138.9±19.0 


2.85±0.06 


25.2±1.2 


0.60±0.07 




g 


69.3±16.8 


3.47±0.11 


2.6±1.7 


0.94±0.21 




ll 


31.7±12.9 


6.22±0.23 


5.1±2.1 


1.21±0.46 




i 


143.6±29.6 


8.42±0.13 


15.5±0.9 


1.36±0.27 




j 


133.8±46.1 


23.67±0.63 


0.1±1.5 


3.66±1.25 


2005.562 


c 


1845.7±53.1 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




(1 


795.9±40.0 


0.58±0.06 


1.9±5.9 


0.26±0.01 




e 


139.4±23.1 


1.65±0.06 


18.8±2.1 


0.70±0.10 




f 


139.9±17.6 


2.91±0.06 


26.5±1.2 


0.67±0.07 




g 


70.3±16.6 


3.57±0.13 


4.1±2.0 


1.18i0.26 




h 


69.0±26.7 


6.43±0.59 


6.3±5.3 


3.09±1.19 




i 


148.5±23.7 


8.55±0.10 


15.1±0.7 


1.32±0.20 




j 


131.1±40.4 


23.84±0.56 


0.6±1.4 


3.71±1.13 


2006.515 


c 


1499.6±30.5 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




(1 


502.0±22.9 


0.74±0.06 


-3.8±4.6 


0.33±0.01 




c 


265.8±24.4 


1.41±0.06 


11.6±2.4 


0.76±0.06 




f 


132.2±15.7 


3.12±0.06 


26.6±1.1 


0.66±0.07 




g 


89.5±23.9 


4.72±0.28 


4.6±3.4 


2.20±0.57 




i 


157.8±23.5 


8.73±0.10 


14.2±0.7 


1.48±0.21 




j 


132.6±40.1 


23.92±0.58 


0.4±1.4 


3.90±1.17 


2007.099 


c 


1378.7±51.5 


O.OOiO.OO 


O.OiO.O 


0.07±0.01 




(1 


553.0±30.8 


1.06±0.06 


-7.3±3.2 


0.35±0.01 




e 


282.7±25.8 


1.56±0.06 


15.3±2.2 


0.56±0.04 




f 


124.7±17.3 


3.26±0.06 


27.2±1.1 


0.70±0.08 




g 


66.8±14.5 


5.15±0.22 


4.8±2.4 


2.09±0.43 




i 


152.5±21.5 


8.82±0.10 


13.7±0.6 


1.48±0.19 




j 


124.4±35.3 


24.04±0.53 


0.6±1.3 


3.77±1.06 


2007.441 


c 


1138.4±60.1 


O.OOiO.OO 


O.OiO.O 


0.07±0.01 




(1 


545.9±20.6 


1.17±0.06 


-8.6±2.9 


0.31±0.01 




c 


310.9±47.4 


1.52±0.06 


15.4±2.3 


0.62±0.09 




f 


109.6±19.5 


3.29±0.06 


27.0±1.0 


0.74±0.11 




g 


67.1±19.4 


5.55±0.35 


5.0±3.6 


2.49±0.70 




i 


139.5±22.7 


8.87±0.12 


13.5±0.8 


1.52±0.23 




j 


121.9±37.0 


24.19±0.65 


0.9±1.5 


4.33±1.30 


2008.333 


c 


2765.9±57.8 


O.OOiO.OO 


O.OiO.O 


0.04±0.01 




(1 


357.8±30.3 


1.61±0.06 


-8.8±2.1 


0.33±0.02 
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Table 9 — continued 



Epoch 


Id. 


Flux 


Core Separation 


P.A. 


Size 






[mjy] 


[mas] 


[degree] 


[mas] 




e 


419.9±35.6 


1.77±0.06 


8.7±1.9 


0.90±0.07 




f 


89.8±15.0 


3.53±0.06 


25.1±1.0 


0.72±0.11 




g 


52.4±12.6 


6.04±0.18 


7.6±1.7 


1.61±0.36 




i 


157.9±27.5 


9.11±0.16 


12.9±1.0 


1.96±0.33 




j 


136.8±32.0 


24.21±0.44 


0.5±1.0 


3.84±0.89 


2008.544 


c 


2928.7±62.2 


O.OOiO.OO 


O.OiO.O 


0.07±0.01 




(1 


225.0±24.5 


1.65±0.06 


-9.3±2.1 


0.26±0.02 




c 


516.2±34.8 


1.74±0.06 


6.3±2.0 


0.99±0.06 




f 


87.5±12.8 


3.56±0.06 


24.9±1.0 


0.76±0.10 




g 


60.6±15.7 


6.01±0.22 


7.8±2.1 


1.84±0.45 




i 


161.2±30.0 


9.15±0.18 


12.9±1.2 


2.07±0.37 




j 


143.0±34.9 


24.32±0.48 


0.6±1.1 


3.96±0.95 


2008.757 


c 


3197.9±39.9 


O.OOiO.OO 


O.OiO.O 


0.05±0.01 




(1 


155.7±12.8 


1.75±0.06 


-9.1±2.0 


0.22±0.01 




c 


456.9±33.9 


1.83±0.06 


6.9±1.9 


0.87±0.06 




f 


83.0±15.9 


3.43±0.09 


27.7±1.4 


1.01±0.17 




g 


42.4±13.6 


5.89±0.15 


4.4±1.4 


1.05±0.29 




i 


160.1±29.7 


8.86±0.16 


13.2±1.0 


1.78±0.31 




j 


140.1±41.5 


24.15±0.57 


0.7±1.3 


3.89±1.14 


2009.153 


c 


2929.0±73.9 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




n' 


842.8±42.8 


0.28±0.06 


8.1±12.1 


0.03±0.01 




(1 


89.2±19.2 


1.80±0.06 


-10.3±1.9 


0.29±0.05 




c 


489.0±57.5 


1.97±0.06 


4.9±1.7 


1.07±0.12 




f 


65.3±25.2 


3.58±0.13 


24.1±2.0 


0.75i0.26 




g 


43.3±16.8 


6.26±0.25 


4.0±2.3 


1.36±0.51 




i 


154.1±48.9 


8.99±0.29 


12.8±1.9 


1.87±0.58 




j 


126.4±53.5 


24.31±0.79 


0.6±1.9 


3.75±1.57 


2009.482 


c 


3125.8±86.9 


O.OOiO.OO 


O.OiO.O 


0.09±0.01 




n' 


166.3±32.2 


0.46±0.06 


2.9±7.4 


0.06±0.01 




(1 


41.9±9.3 


1.78±0.06 


-9.4il.9 


0.11i0.02 




c 


478.4±48.5 


2.11±0.06 


5.0±1.6 


1.24±0.12 




f 


39.4±9.5 


3.65±0.07 


24.7±1.1 


0.67±0.14 




g 


38.4±11.5 


6.47±0.12 


6.1±1.0 


0.89±0.23 




i 


166.0±38.6 


8.83±0.26 


12.9±1.7 


2.32±0.52 




j 


134.3±40.5 


24.27±0.57 


1.0±1.3 


3.84±1.14 
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